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THE ELECTRON BOLTZMANN EQUATION IN A PLASMA 


GENERATED BY FISSION PRODUCTS 
by 

H. A. Hassan and Jerry E. Deese 
North Carolina State University 
Raleigh, North Carolina 

SUMMARY 

A Boltzmann equation formulation is presented for the determination of 
the electron distribution function in a plasma generated by fission fragments. 
The formulation takes into consideration ambipolar diffusion, elastic and in- 
elastic collisions, recombination and ionization and allows for the fact that 
the primary electrons are not monoenergetic . Calculations for He in a tube 
coated with fissionable material show that, over a wide pressure and neutron 
flux range, the distribution function is non-Maxwe Ilian but the electrons are 
essentially thermal. Moreover, about a third of the energy of the primary 
electrons is transferred into the inelastic levels of He. This fraction of 
energy transfer is almost independent of pressure and neutron flux but in- 
creases sharply in the presence of a sustainer electric field. 

INTRODUCTION 

1-3 

Recent experiments in CO, Xe-He and Ne-N^ mixtures have demonstrated 
direct nuclear pumping. At a time when the quest for high power lasers is at 
its peak, such a development is rather significant because it affords the 
means of uniform pumping of large volume high pressure gases. 

Although the inversion mechanisms were not discussed in any detail in 
Refs. 1-3, it is generally believed that both the high energy heavy particles 
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and the low energy electrons play a role in population inversion. ’ An es- 
timate of the extent of the role played by the heavy particles and electrons 
in nuclear pumping requires, as a first step, the determination of their re- 
spective energy distributions. Using a monoenergetic source of high energy 
primary electrons, Wang and Mlley^ calculated the electron distribution func- 
tion using Monte Carlo techniques. Later, Lo and Miley^ employed a simplified 
version of the Boltzmann equation, and the results compared favorably with the 

Monte Carlo calculation. In their study of electric discharges in gases, 
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Thomas and Thomas and Thomas have shown that predictions based on the nu- 
merical solution of the appropriate Boltzmann equations are in good agreement 
with those obtained using Monte Carlo methods. Moreover, computations using 
the Boltzmann equation formulation were carried out at a considerable saving 
in computer time. Because of this and because of the role played by the 
electrons in nuclear pumped lasers, a Boltzmann equation formulation is pre- 
sented for the calculation of the electron distribution function. In this 
formulation, the effects of ambipolar diffusion, elastic and inelastic colli- 
sions, two and three body recombination and secondary ionization are taken 

into consideration. Moreover, because the primary electrons generated by the 

9 

fission fragments are not monoenergetic , the present formulation allows for 
a source of primary electrons whose distribution is calculated by a procedure 
similar to that of Ref. 9. 

The resulting nonlinear differential-difference-integral equation is 
solved for a He plasma generated by fission fragments in the presence and ab- 
sence of externally applied electric fields over a wide range of pressure and 
neutron flux. The results show that the electrons are essentially thermal, 
but the distribution function is far from Maxwellian. 
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ANALYTICAL FORMULATION 


The experiments of Refs. 1-3 employed tubes coated with fissionable ma- 
terial. Under neutron bombardment, fission fragments emerge from the coating 
and enter the gas- The ensuing energy transfer results in ionization and ex- 
citation of the background gas. For the high pressures of interest, the 
plasmas generated by the fission fragments are slightly ionized.^ Therefore, 
a Lorentz gas approximation will be employed. 

Letting the electron distribution function be expressed as: 


v.f . „ 

f = f^(v,x^,t)H V = i = 1,2,3 


( 1 ) 


where is the velocity component, then the governing equations for f^ and 
f^ are given by 
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and e is the electronic charge, m is the electronic mass, M is the mass of 
the heavy particle, T is the gas temperature, N is the gas number density, 
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is the momentum transfer cross section, is the excitation energy, 

is the excitation cross section and (8f /8t) is the source term resulting 

o c 

from secondary ionization, recombination and production of electrons by fis- 
sion fragments. 

Steady state conditions are to be considered. Equation (3) can be solved 
then directly for f^ 
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Substituting equation (6) into equation (2) one finds 
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In general E^ has values along and normal to the axis. The component 

along the axis, E^, is that due to the applied electric field and is usually 

given. The component normal to the axis, E^, is obtained from consideration 

of the diffusion process in the tube. 

Before one can attempt the solution of equation (7), one needs to specify 

E and (3f /3t) . In the presence of ambipolar diffusion E is determined 
z o c z 

from the requirement that the ion flux is equal to the electron flux in the 

normal direction. The ion or electron flux, F. , is defined as 

X, s 
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where s refers to the ions or electrons. 


Using equation ( 6 ) with e replaced 


by -e , one obtains 
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where D is the diffusion coefficient and y is the mobility. The mobility is, 
by definition, a positive quantity; thus. 
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The equation governing n^ follows from equation (7) by multiplying by 
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41TV and Integrating from 0 to Infinity. Carrying out the Indicated opera- 
tion, one obtains 
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To carry out the Integration In equation (14), one needs to determine 

(9f /9t) . This term consists of a secondary ionization term, a production 
o c 

term resulting from the fission fragments and a recombination term. These 
processes are rather complex, making it necessary to employ crude approxima- 
tions for their representation. Thus, the procedure employed by Chan and 
12 

Moody for representing secondary ionization will be employed. Letting 
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the contribution to (9F^/9t)^ resulting from secondary ionization can be 
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written as 


where 
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is the ionization potential and is the ionization cross section. Using 
(17) in (16), one finds 
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To estimate the production term resulting from the fission fragments, a 
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procedure similar to that employed by Guyot, Mlley and Verdeyen will be 
used. The fission fragments are assumed to fall into two groups: a light 
group with an average mass number of 96, an average initial charge of 20e and 
an average energy of 98 MEV; and a heavy group with an average mass number of 
140, an average initial charge of 22e and an average energy of 67 MEV. 

Consider the case of a tube coated with fissionable material and bom- 
barded with neutrons. For the case where the range of the particle in the 
coating is less or equal to the thickness of the coating, the number of heavy 
particles k born with energy E , and having an energy E at location z per 

O y Ki 
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unit surface per unit time and per unit energy is chosen as 
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where 


g.k 
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o,k 
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X is the range of the particle in the gas, T is the thickness of the coating 
S 

and S is the source rate of charged particles per unit volume and is given by 


S = Og (j)/2 


where (j) is the neutron flux, is the number density of the target material 
and is the fission cross section. For fission fragments n - 1/2. 

For a slab model any quantity Q(y) at a distance y from one of the planes 
is given by 
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The contribution from two plane parallel sources at a distance 2a apart gives 


= Q(y) + Q(2a - y). 


( 22 ) 


The number of electrons created per unit volume, per unit time, and per 
unit energy with a kinetic energy between E and E + dE at location z by the 
k*"^ heavy particle is given by 
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and N Is the number of electrons in the s shell of the gas atom, e. is 

S to > 

the ionization potential, Z ,= 20 or 22 and m == 1 for fission fragments. 

O ^ K. 

The expression for follows from Ref. 13. The desired contribution to 
(3F^/9t)^ is obtained from equations (22), (23) and (24) as 


g(y>£) = 1 {f,,(y.e) + fv^^a - y,e)}. 
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Two and three body recombinations are considered. For a recombination 
reaction of the type 

+ e ^ A 


the contribution to (9F /9t) can be written as 
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three body recombination, the contribution to (9F /9t) will be chosen as 
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where a is the three body recombination coefficient. The reasoning that led 
to equation (27) is similar to that employed in deriving equation (16). The 
choices indicated in equations (16), (26) and (27) satisfy both particle and 
energy conservation. 

f 2 

Using equations (16), (25), (26) and (27), the term 47T v (8f /dt) dv 
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Because g(z,e) does not depend on a simple way on position and because, for 
the pressure and dimensions under consideration, it is well approximated by 
its value at the axis, it is assumed when solving equation (14) that 

3(z) = 3(0). (32) 

Although the above assumption simplifies the calculations, it is not required 
for carrying out the solution. The boundary conditions for equation (14) can 
be expressed as 

dn 

-T-^ = 0 at z = 0 (33) 
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and, at the wall, a Schottky boundary condition is employed, i.e. 

n = 0 at z = a. (34) 

e 

With n determined from equation (14), E follows from equation (12). 
e ^ z 

METHOD OF SOLUTION 


The function f (v,z) = H (e,z) 
o o 
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where the n.’s represent an orthogonal set. This set may be obtained from a 

suitable linear combination of functions ^ which are 

linearly independent for arbitrary r. Ref. 14. By definition, the electron 
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where n (z) is determined from equation (14). The governing equation for h, 
e ^ 

is determined from equation (7) by substituting first for E^ and (9f^/9t)^ and 
then multiplying by n (z) and integrating from zero to a. The resulting 
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equation for hj^ depends on ... This means that before one can deter- 

mine h, , one has to determine h , h, . . . h, , . 
k o 1 k-1 

In this work, it is assumed that the solution may be approximated by the 
first term. Using the above procedure, one finds that h^ is given by the 
following equation. 
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As may be seen from equations (14) and (30), the O's are not independent. 
Thus, multiplying equation (14) by n^ and integrating between o and a, one 


finds 
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Equation (42) follows also from equation (39) by Integrating with respect to 
E from 0 to Infinity. 
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it is seen that equation (39) can be expressed as two first order equations 
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Equations (44) and (45) can be integrated using a Runge-Kutta or a predictor- 
corrector method starting at some £ = £ ^^ ^. Before one can start the integra- 
tion one has to determine a,; these quantities depend on n (z) which is gov- 

J ® 

erned by equations (14) and (30). Because some of the quantities appearing 
in equation (30), i.e. D. and a, depend on h , while others, i.e. S„ and a , 
depend on T^, the electron temperature which, in turn, depends on h^, the 
problem under consideration is nonlinear and an iterative procedure is re- 
quired to determine the solution. 

As a result of equation (42) 

G = 0 at £ = e (47) 

max 


thus, to start the integration, one needs to assume a, , y , D and the 

K, O G 6 

value of h at £ = £ . For the assumed a, y , D , and A„ , the value of 

o max e e K,o 

h at £ = £ is determined from the requirement indicated by equation (37) . 
o max 

Using the calculated h^(c) new values of D^, y^, A^ ^ and a are calculated 
from equations (10), (11), (29) and (31) and are used to recalculate . The 
above procedure is repeated until convergence is achieved. 
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As is seen from the above, the calculation of h^ requires a rather 
lengthy iterative procedure. Unfortunately, the convergence of the above 
method becomes increasingly difficult as the electric field decreases and/or 
as the pressure Increases. This is because the coefficients of the highest 
derivatives in equation (39) become rather small thus necessitating extremely 
small integration steps. Because of this, the method of composite expansions 
was used in integrating equation (39) when was set equal to zero. 

In this method the function h (u) is written as 
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h^ = h^(n;(S) = 'J^(n;<S) + x(Uq;5) 
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where 6 is the largest of the two parameters (2m/M) and (ey/e^)!^ and 
£ 

= ri/6 is the inner variable with £ > 0 being determined according to the 
procedure outlined below. The function X(U^>5) is negligible outside the in- 
ner region, i.e. when ri^ ->■ 

The equations for are determined by setting 

h = Z 6^ ip (49) 
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in equation (39) and equating equal powers of 6. When this is done is de- 
termined from all terms in equation (39) that are not multipled by 6, i.e. 
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The equations for X^(O^) are determined by first transforming equation (39) 

from ri to ri^ and choosing Si so that the coefficient of the highest derivative 

does not depend on 6 . This procedure was not used here; instead, equations 

(44) and (45) were integrated starting at r| = 0 using the condition G(0) = 0 

and an assumed value of h^(0). The quantity h^(0) was determined from the 
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where 6 and 6 are preselected small numbers, 
m n 

As a check on the accuracy of the Integration procedure and in order to 
determine the manner in which the primary electrons dispose of their energy, 
an electron energy equation is derived. This equation is obtained by multi- 
plying equation (30) by f| and integrating from zero to infinity. The result- 
ing equation can be written as 
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The above equation shows that the energy received by the electrons from fis- 
sion fragments, electric field and three body recombination is equal to, re- 
spectively, the energy lost by diffusion, elastic collisions, inelastic col- 
lisions, two body recombination and secondary ionization. 

RESULTS AND DISCUSSION 

The solution of equation (39) is carried out for a He plasma. Before a 

solution can be obtained one needs to know the momentum transfer, excitation 

and ionization cross sections of He. The momentum transfer cross section up 

1 6 

to an electron energy of 6 ev is taken from Crompton et al. ; above 6 ev it 

9 

is assumed that the collision frequency V is a constant equal to 2.4 x 10 

sec \ Ref. 12. All helium excited states with a principal quantum number of 

5 or less are included in the calculations and the excitation cross sections 

for these states are taken from Refs. 17-25. The ionization cross section is 

26 

taken from Rapp and Englander-Golden. 

At the high pressures of interest, the recombination process in noble 
gases is complicated by the formation of molecular ions. The theory presented 
here allows for radiative and three body recombination. However, at high 
pressures the reactions 

He + 2 He ^ He^ + He (53) 

e + He^ + He -> He + 2 He (54) 

where * designates an excited state play a dominant role since the forward 

rate for the first reaction^^ is about 10 cm^/sec at 300°K while the for- 

28 “27 6 o 

ward rate for the second reaction is about 2 x 10 cm /sec at 300 K. Be- 
cause of the rapid conversion of He to He 2 at high pressure it is seen from 
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the second reaction that an effective two body recombination coefficient can 
be defined whose approximate value is 

B ~ 2 X 10 n cm^/sec. (55) 

R ri6 

29 

For reactions of the type indicated in equation (54) , Biondi suggested a 
forward rate coefficient of the form 

= 10"^^^”^^(Te/300)“^'^^ cm^/sec. (56) 

Because of the above uncertainties, the calculations presented here assume a 
-9 

value of of 10 at a pressure of 100 Torr; a value which was employed in 
R 

Ref. 5. At .’other pressures is scaled according to the estimate indicated 

R 

in equation (55) . 

2 

The ion mobility is taken as 10.4 cm /V sec, which is appropriate for a 
27 30 

He ion. ’ Because the ion mobility does not play a significant role in 
the present calculations, no attempt was made to allow for the presence of 
other He ions. 

The calculations were carried out for a coating of U with the neutron 

3 o 

flux ranging from 3.8 x 10^^ to 7.6 x 10^^ neutrons/cm^ sec at a temperature 
of 300°K and pressures of 100 and 760 Torr. The range of parameters was cho- 
sen to ensure that the Lorentz gas approximation employed here is not violated. 
The spacing between the walls, i.e., 2a, is taken as 3.7 cm. 

A typical plot of the distribution function of the primary electrons, 
g(e), is shown in Fig. 1. The plot demonstrates that the energy spectrum of 
the primary electrons is far from monoenergetic . Figure 2 shows the effect 
of the neutron flux on the electron distribution function F^ (See equation 
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(15)). Examination of equation (39) shows that at high energies 


a: g(n) “ g(n) . 


( 57 ) 


The distribution function of the primary electrons, g, is proportional to the 
neutron flux. On the other hand, conservation of particles, equation (42), 
is such that is approximately inversely proportional to the square root of 
the neutron flux; thus, F^ Increases with the neutron flux at high energies. 
In the low energy region, the distribution function approaches a Maxwellian 
at the gas temperature. 

The effect of pressure is discussed next. As a result of increased col- 
lisions, the high energy particles are depleted at a faster rate thus result- 
ing in a reduced distribution at higher pressures. Because F^ vs e is pre- 
sented on a log-log scale and because F is normalized, calculations of F at 
100 and 760 Torr show that the effect is small and, therefore, does not war- 
rant a separate plot. 

It is evident from Fig. 2 that the distribution function is not Maxwel- 
lian. However, calculation of the electron temperature, i.e. 


T 

e 


_ 2 _ 

3k 


e F de, 

o 


(58) 


shows that the electrons are essentially thermal. Assuming a Maxwellian dis- 
tribution function at this temperature would result in a substantial reduction 
of both electron excitation and ionization. This suggests that assuming the 
distribution function to be Maxwellian at some temperature different from the 
gas temperature is not a good approximation. Ref. 4. 

The effect of an externally applied electric field is shown in Fig. 3. 

It is seen from the plot that the distribution function approaches a Druyvestyn 
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distribution at low energies. At high energies the distribution is similar 
to that for a zero electric field. 

Figure A shows the manner in which energy is transferred from the primary 
electrons. Over a fourth of the energy is transferred to the excited states 
and this fraction is essentially independent of the range of neutron flux and 
pressure considered here. Energy transfer from electrons in the presence of 
an electric field was also investigated. The results indicate that over 90% 
of the energy is transferred to the He excited states. Unfortunately, for the 
E/P considered, namely .1 V/cm Torr, the energy received by the electrons from 
the electric field far outweighed the energy received from the fission frag- 
ments. In spite of this, the above clearly illustrates the beneficial effect 
of employing a sustalner electric field in conjunction with nuclear induced 
plasmas for the generation of high power lasers. 

Because the analyses of Refs. 5 and 6 employ monoenergetic primary elec- 
trons and do not employ the experimentally measured cross sections employed 
here, direct comparison with the present work is not possible. Our results 
are, however, in qualitative agreement with their results. Based on present 
and some earlier unpublished calculations, we feel that the usual characterl- 
zation of diffusion , namely, ignoring the ambipolar part and using a cosine 
(or Bessel function) representation to approximate the spatial derivative of 
the distribution function, is not adequate for fission generated plasma (See 
equations (14) and (30)). 

In this work, the solution of the Boltzmann equation is approximated by 
the first term at the series solution indicated in equation (35) . This ap- 
proximation implies that the electron temperature is constant throughout and 
thus is consistent with the assumption of constant gas temperature. This 
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explains the reason why energy losses by conduction are not included in Fig. 

3. Therefore, as long as the assumption of constant temperature is adequate, 
which is the case at high pressures, there is no need to carry the computa- 
tions beyond the first term of the assumed series. 

CONCLUDING REMARKS 

The formulation presented here for the calculation of the electron dis- 
tribution function in plasmas generated by fission fragments is quite general 
and may be used for any gas or gas mixture. Results based on this formula- 
tion for He show that a large fraction of the energy of the primary electrons 
is transferred into the excited states. Moreover, when such plasmas are sub- 
jected to a sustainer electric field, sharp increase in the excitation rates 
results. Therefore, it appears that the major contribution of fission frag- 
ments is to provide the means for generating plasmas at high pressures. These 
plasmas can, in turn, be employed in conjunction with sustainer electric 
fields to generate high power lasers. 
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Figure 4. Fractional power transferred from the 
high energy electrons. 
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